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ELECTRICAL  PROPERTIES  OF  CONDUCTING  GALLIUM  ARSENIDE 


1.1  Introduction 

Tho  two  basic  parameters  determined  in  electrical-measurement  studies 

are  l  lie  carrier  concent  rat  ion  (n  or  p)  and  the  carrier  mobility  (p  or  u  ). 

n  p 

(Hereafter,  we  will  refer  only  to  n  and  unless  we  are  specifically  dealing 
with  p-tvpe  samples.)  The  most  common  experiment  used  to  ascertain  n  and  u 

n 

is  the  Hall-effect,  along  with  conductivity.  To  use  the  Hall-effect  in  the 
most  accurate  manner,  it  is  necessary  to  know  the  Hall  r-factor,  defined  by 
n  =  r/eR,  where  R  is  the  Hall  coefficient.^  Then,  once  n  is  known  accur¬ 
ately,  it  must  be  properly  related  to  the  real  quantities  of  interest, 
namely,  the  donor  and  acceptor  concentrations  (N  and  N  ,  respectively), 

V)  A 

and  their  energies  with  respect  to  the  relevant  band  (E  and  E  ,  respec- 

L)  A 

tivoly).  However,  at  the  beginning  of  this  contract  it  was  clear  that 
neither  of  these  subjects  was  being  handled  well  in  common  practice.  For 
example,  nearly  all  researchers  simply  dealt  with  the  first  problem  by 
assuming  r  =  1,  and  with  the  second,  by  assuming  K  =  1,  4,  or  -i  in  the 
distribution  function  for  donors  or  acceptors: 

f(S)  =  1/  [1  +  K  exp  (&-ii»  )/  kT]  (1-1) 

r 

where  is  the  energ-1'  of  an  electron  and  ^  is  the  Fermi  energy.  However, 

r 

even  the  form  of  Eq.  1-1  is  wrong  if  a  bound  center  has  more  than  one  charge 
state  within  the  band  gap,  or  important  excited  states.  Thus,  we  realized 
that  some  theoretical  work  was  needed  before  any  serious  conclusions  could 
be  drawn  from  our  electrical  measurements.  Also,  hole  transport  theory  is 
even  in  much  worse  shape  than  electron  transport  theory,  and  we  were  able 
to  make  some  improvements  here,  although  nothing  fundamental. 

We  applied  these  concepts  to  a  number  of  problems  in  conducting 
Ga  As ,  i  in  I  iid  iiij; :  (I)  a  detailed  analysis  of  low-compensation  vapor-phase- 

epilaxial  (Vi’E)  GnAs;  (2)  the  identification  of  two  prominent  centers  at 


1 


H  — 0 .IS  cV,  and  K  — i) .43  eV,  as  pure  defects;  and  (1)  the  development  of 

1*  l* 

mobility  profiling  techniques  for  conductive  FET  layers.  These  problems 

will  he  discussed  in  detail  below. 

1.2  Statistics  of  Multi-Charge  Centers  in  Semiconductors 

Most  researchers  seriously  involved  with  electrical  measurements, 

and  in  particular  with  temperature-dependent  Hall-effect  measurements, 

have  experienced  trouble  with  choosing  or  determining  a  degeneracy  factor, 

i.e.,  the  "K"  of  Eq.  1-1.  Thus,  we  dealt  with  this  problem  and  derived 

2 

the  tol lowing  general  formulas: 


k  v1  m  = 


l  +  >: 


(1-2) 


Vm  '  \iv  >nK 


2. '  ,m'  I2,m 


n  -  +  £  <«k-Onktm 

k,  2,m 


(1-3) 


where  is  the  concentration  of  type-k  donors  or  acceptors,  2^  is  the 

number  of  ionizable  electrons  (with  less  than  band-gap  energy)  on  a  center 

of  type-k,  "m"  denotes  excited  states,  is  the  degeneracy  of  the  k2m^— ■ 

state,  n,  a  ~  n  /kT,  and  1  =  0,  1,  —  2,  .  The  restriction  on  the  summa- 
kv.m  kv.m  k 

tion  means  that  2+2,  and  mtm  ,  a_t  the  same  time.  Equations  1-2  and  1-3 
are  completely  general  for  any  number  of  donors  and  acceptors  with  any 
numbers  of  ionizable  electrons  and  excited  states.  Note  that  Eq.  1-3  is 
simply  the  charge-balance  equation.  Many  examples  of  the  application  of 
these  equations  are  given  in  References  1  and  2. 


V  •,  V .  jv'j 


I.)  K1  citron  Transport  Theory 

To  properly  use  liquations  1-2  and  1-3,  wo  must  know  n  (or  p,  for 
p-tvpo  samplos),  as  a  function  of  Lompo rat ure.  From  a  Hall-effect  measure¬ 
ment  we  get  n  =  r/eR,  where  R  is  the  measured  quantity,  but  r  is  unknown. 

Most  researchers  simply  set  r  =  1,  which  is  justified  for  rough  measure- 

<  < 

ments,  since  1  -  r  ~  1.3  in  the  usual  GaAs  samples.  However,  to  do  better, 

2 

we  must  either  measure  R  at  very  high  magnetic  fields  (such  that  U(cm  /V-sec) 
£ 

B  (gauss)  10  ) ,  or  else  calculate  r  from  mobility  measurements,  fitted 

to  scattering  theory.  For  relatively  pure  GaAs,  in  the  room-temperature 

range,  polar  optical-mode  scattering  is  dominant,  and  it  is  necessary  to 

solve  the  Boltzmam  transport  equation  by  numerical  means.  We  have  adopted 

3  4 

a  program  written  by  Nag  ,  who  uses  Rode's  iterative  solution  of  the 
Boltzmann  equation.  Modifications  to  Nag's  program  include  the  considera¬ 
tion  of  neutral  impurities  in  both  scattering  and  screening.  The  use  of 
this  program  to  fit  data  is  presented  in  a  later  section.  Recently,  we 
have  also  modified  the  ionized-impurity  scattering,  by  incorporating  the 
work  of  Meyer  and  Bartoli'*.  Various  unknown  parameters,  such  as  the 
acoustic-deformat  ion-potential  constant,  have  been  determined  by  fitting 
to  two  pure  VPF  layers,  grown  at  AFWAL/AADR  by  Peter  C.  Colter.  Thus,  we 
feel  that  our  mobility  fitting  program  is  probably  one  of  the  most  complete 


in  existence,  at  the  present  time.  Much  of  the  computer  programming  for 
this  project  has  been  carried  out  by  John  R.  Sizelove,  AFWAL/ADO. 


1.4  Hole  Transport  Theory 

The  following  abstract  is  reprinted  from  Ref.  6. 

We  have  used  a  two-band  model  (heavy  and  light  holes)  to  calculate 
the  transport  properties  of  p-type  (JaAs.  The  scattering  mechanisms 
included  are  acoustic-mode  deformation  potential,  acoustic-mode  piezo¬ 
electric  potential,  polar-  and  nonpolar-mode  deformation  potential,  ionized 
impurity,  and  space  charge.  Interband  scattering  is  included  explicity  for 
the  optical  phonons  and  phenomenologically  for  the  acoustic  phonons.  The 
intraband  polar  optical-mode  scattering,  for  which  a  relaxation  time  cannot 
be  defined,  was  calculated  by  using  the  numerical  method  of  Fletcher  and 
Butcher.  The  acoustic  deformation-potential  parameter  and  the  coupling 
coefficient  for  interband  scattering  were  calculated  by  fitting  the  theory 
to  Hall-mobility  data  for  both  pure  and  doped  samples.  We  have  determined 

ion ized-impuritv  and  space-charge  contributions  for  two  of  our  samples,  one 

4+  3+ 

doped  with  Cr  and  exhibiting  an  0.33-eV  activation  energy  (Cr  /Cr  ),  and 
the  other  heat-treated  and  exhibiting  an  0.14-eV  activation  energy.  The 
transport  fit  was  carried  out  self-consistently  with  a  statistical  fit  of 
the  p  vs  T  data,  to  determine  the  relevant  donor  and  acceptor  concentrations 
Individual  plots  of  the  mobility  contributions  from  the  various  scattering 
mechanisms  are  also  presented.  Finally,  useful  plots  of  the  Hall  mobility 
and  Hall  r  factor,  as  functions  of  carrier  density  and  compensation  ratio, 
are  provided  to  aid  others  in  the  analysis  of  their  samples.  An  important 
conclusion  of  this  study  is  that  the  light-hole  band  strongly  influences  the 
transport  properties,  even  though  its  density  of  states  is  small  compared  to 
that  of  the  heavv-hole  band. 


1.5  Electrical  Properties  of  Low-Compensation  GaAs 
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It  was  found  by  P.  C.  Colter,  of  A  is  laboratory,  that  very  low 
receptor  concentrations  could  be  obtained  in  VPE  GaAs  by  baking  the  Ga 
source  prior  to  growth  to  reduce  Zn  contamination,  and  by  taking  other 
precautions  to  reduce  deep  levels.  Because  of  the  uniqueness  of  these 
samples,  a  very  careful  electrical-measurement  study  was  in  order.  First 
we  measured  R  vs.  T  in  an  automated  Hall-effect  system,  developed  earlier^. 
Then  we  determined  n  from  the  relationship  n  =  r/eR,  where  r(T)  was 
calculated  from  our  mobility-fitting  program.  Finally,  we  fitted  n  vs.  T 
to  the  following  charge-balance  equation,  developed  from  Equations  1-1  and 


n  +  N  =  Nn  /  (1  +  n/(f>) 
A  o 


where 


(1-4) 


d>  =  2(27Tm  *kT) 
n 


o/k  3/2  -  E  n/T 


(1-5) 


Here  the  donor  energy  is  defined  by  =  E^  -  aT,  and  gQ  and  g^are  the 
degeneracies  of  the  unoccupied  and  occupied  donor  states,  respectively. 
From  the  n  vs.  T  fit  we  get  the  parameters  Np,  N^,  E^,  and  exp 

(a/k) .  We  also  get  N  from  the  mobility  fit.  The  experimental  and 

A 

theoretical  curves  of  y  and  n  are  presented  in  Figures  1-1  and  1-2, 

n 

respectively,  and  the  fitted  parameters,  in  Table  1—1.  Both  the  acceptor 
concentration,  and  the  compensation  ratio  (N^/NQ) ,  are  the  lowest  ever 
reported  for  GaAs. 


For  further  details  of  this  work,  see  Ref.  8. 
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l.b  Defect  Centers  Ln  CaAs 


As  crystal  growers  have  been  able  to  reduce  impurity  contamination 
in  GaAs,  the  intrinsic  defects  have  become  more  important.  Although  many 
such  defect  energy  levels  have  been  found  in  electron-  or  ion-irradiated 
samples,  none  of  these  has  ever  been  proven  dominant  in  as-grown  material. 
(There  is,  however,  now  a  strong  belief  that  the  dominant  deep  donor 
(EL2)  in  undoped,  semi-insulating  LEC  material,  is  indeed  the  anion  antisite 
defect  As  .)  In  general,  n-type,  bulk,  as-grown  CaAs  is  dominated  by 
one  of  three  levels:  ( L)  the  shallow,  hydrogenic  levels,  such  as  Si,,  or 
Se  ;  (2)  an  0.15  eV  level;  and  (3)  an  0.43  eV  level.  We  have  shown  by  a 
combination  of  electrical  and  analytical  techniques  that  the  latter  two  are 
pure  detects,  not  associated  with  any  impurities.  The  technological 
importance  of  this  result  is  that  further  reductions  of  donor  and  acceptor 
concentrations  in  bulk  GaAs  will  have  to  be  concerned  with  defects  as  much 
as  impurities. 

Further  details  may  be  found  in  References  9  and  10. 

1.7  Mobility  Profiling  in  FET  Layers 

Although  carrier-concentration  profiles  in  semiconductor  layers 

have  been  routinely  extracted  for  years,  by  capacitance-voltage  (C-V) 

measurements  on  Schottky  barriers,  it  is  much  more  difficult  to  measure 

mobility  profiles.  Two  methods  which  have  been  developed  include:  (1) 

the  Hall-effect  method,  in  which  a  reverse-biased  Schottky  barrier  is 

used  to  successively  deplete  the  conductive  layer^S  and  (2)  the  FET 

geometric  magnetoresistance  (GMR)  method,  in  which  the  magnetic-field 

dependence  of  the  source-drain  (SD)  current  is  measured  as  the  reverse- 
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bias  gate  voltage  is  increased,  in  steps  ’  .  We  have  developed  a 

fully  automated  apparatus  to  carry  out  both  of  these  techniques.  The 


UMR  technique  is  simpler  in  that  only  one  parameter  must  be  measured, 
namely  the  S-D  current,  as  long  as  gate  current  is  negligible.  However, 
to  obtain  data  close  to  the  surface,  it  is  necessary  to  forward-bias  the 
FET,  and  then  gate  current  can  become  important.  We  have  developed  an 
FET  circuit  model  in  which  gate  current  is  properly  taken  into  account. 
If  it  is  impossible  to  measure  the  gate  current,  then  there  are  ways  to 
minimize  its  effect,  e.g.,  by  averaging  1^  readings  for  and  ~VgD* 

In  Fig.  1-3  we  show  the  carrier  concentration  and  mobility  profiles  in 
a  FAT  FET,  for  which  the  gate  current  is  actually  larger  than  the  S-D 
current,  in  the  most  forward-bias  mode.  It  is  seen  that  the  mobility 
correction  is  very  large  in  this  region.  With  these  corrections,  we  are 
now  able  to  get  accurate  mobility  data  much  closer  to  the  surface  than 


before. 


Fig.  1-3  The  electron  concentration  and  mobility  for  a  semi- 

insulating  GaAs  substrate  implanted  with  4  x  10  ^/cm  , 
100  keV  Si  ions,  and  annealed  at  850°C  for  15  min. 


2.0  ELECTRICAL  PROPERTIES  OF  SEMI-INSULATING  GALLIUM  ARSENIDE 


2.1  Introduction 

A  big  advantage  in  the  development  of  GaAs  IC  technology  is  the 

existence  of  semi-insulating  (SI)  forms  of  GaAs,  namely  those  containing  Cr, 

or  the  defect,  EL2.  Thus,  it  is  possible  to  grow  good  single-crystal 

epitaxial  layers,  or  form  a  conductive  layer  by  implantation,  and  not  have 

to  worry  about  the  "platform,"  namely  the  S  I  GaAs.  Unfortunately,  this 

optomistic  picture  has  not  been  fully  realized,  as  yet,  for  various  reasons. 

The  importance  of  the  S  I  GaAs  problem  is  demonstrated  by  the  calling  of 

international  conferences  on  the  subject  in  1980,  1982,  and  in  1984  (coming 

up).  Thus,  we  have  devoted  a  large  part  of  our  effort  to  this  subject  and 

have  made  both  theoretical  and  experimental  contributions.  The  highlights 

3+  4+ 

include:  (1)  discovery  of  a  donor  transition  (Cr  /Cr  )  in  Cr-doped  GaAs; 

(2)  showing  that  an  0.1  eV  "conversion"  center  is  not  always  Mn,  as  commonly 
believed;  and  (3)  correlating  SI  substrate  properties  with  the  device  (FET) 
properties. 

2.2  Magneto-Hall  and  Magnetoresistance  Coefficients  in  Semicon¬ 
ductors  with  Mixed  Conductivity 

The  electrical  properties  of  Cr-doped  GaAs  are  complicated  by  the 

existence  of  both  hole  and  electron  conductivity,  since  the  Fermi  level  is 

near  the  middle  of  the  band  gap.  Thus,  four  unknowns  exist:  n,  p,  y^, 

and  Up.  The  problem  can  be  solved  exactly  by  measuring  the  second-order 

magnetic  field  dependences  of  the  conductivity  and  Hall  coefficient,  if_ 

all  of  the  magnetic-field  dependence  of  these  parameters  is  due  to  mixed- 
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conductivity  effects.  However,  if  single-carrier  magnetic-field 

dependences  are  also  large,  then  the  problem  is  much  more  complicated. 

2  2 

We  have  developed  detailed  expressions,  to  order  y  B  ,  for  the 


magnetoresistanco  and  magneto-Hall  coefficients,  allowing  both  single-carrier 
and  mixed-carrier  effects.  The  expressions  are  too  lengthy  to  be  presented 
here,  but  further  details  may  be  found  in  Ref.  15. 

4+  3-f- 

2.3  Discovery  of  the  Cr  /Cr  Donor  Transition  in  GaAs 

Temperature-dependent  Hall-effect  measurements  on  two  Cr-doped  GaAs 

-4 

samples  show  a  dominant  center  at  =  0.324-1.4  x  10  T  eV,  with  respect  to 
the  valence-band  edge.  By  comparison  with  secondary-ion  mass  spectroscopy 
measurements  of  the  Cr  concentration,  and  recent  EPR  measurements  of  the 

2+  3+  i | 

Cr  ,  Cr  ,  and  Cr  concentration  in  several  samples,  it  is  shown  unambigu- 

4+  3+ 

ously  that  this  energy  describes  the  Cr  ~*Cr  transition.  This  is  the  first 

4+ 

conclusive  evidence  for  a  charge-state  transition  involving  Cr  in  GaAs. 

The  technological  importance  of  this  result  is  that  to  produce  p-type  material 
say  by  implantation  of  acceptors,  in  Cr-doped  material,  we  must  have 
Na  >  ND  +  NCr,  not  just  Na  >  N^,  which  would  be  the  case  if  Cr  did  not  have 
a  donor  transition. 

Further  information  may  be  found  in  Ref.  16. 

2.4  The  0.1  eV  Conversion  Center  in  GaAs 


The  existence  of  a  serai-insulating  (SI)  form  of  GaAs  has  been 
essential  for  the  development  of  much  of  the  present-day  GaAs  device  tech¬ 
nology."^  A  well-known  problem  with  SI  GaAs,  however,  is  the  "conversion" 
phenomenon,  in  which  a  thin,  conducting  surface  layer  is  formed  upon  heating 
to  the  temperatures  required  for  expitaxial  growth,  or  ion-implantation 
annealing.  One  of  the  most  common  types  of  conversion  involves  a  p-type 
layer  which  exhibits  a  thermal  (Hall-effect)  activation  energy  of  0.1  eV, 
and  a  photoluminescence  line  at  1.41  eV.  Such  a  layer  is  formed,  e.g.,  by 
annealing  SI  GaAs  in  H„  gas  at  700-800°C.  The  relevant  acceptor  has  been 


J  / 


attributed  to  Mn  ,  ' 8-2  *  SI  -  V  ,  24-27  C  -  V  ,  2(’"28  and  V  .29  In 
an  attempt  to  study  this  problem  further,  Yu  and  Park2^  have  compared  the 
shapes  and  temperature  dependences  of  the  respective  1.41-eV  emission  lines 
from  Mn-implanted  samples,  and  unimplanted,  annealed  SI  samples.  The 
results  were  nearly  identical.  Thus,  we  feel  that  Mn  represents  the  most 
likely  cause  of  the  1.41-eV  emission  in  the  majority  of  annealed  samples 
investigated  to  date,  although  the  other  possibilities,  mentioned  above. 


cannot  be  ruled  out  in  all  cases. 


We  report  here  a  rather  extensive  investigation  of  two  semi-insulat- 
ing  liquid-encapsulated  Czochralski  crystals,  one  undoped  (sample  F)  and  the 
other  Cr  doped  (sample  G) .  Both  samples  were  annealed  at  750°C,  face  up, 
for  15  min.  in  an  H»  ambient,  to  simulate  vapor-phase  epitaxial  (VPE)  growth 
conditions.  The  surfaces  appeared  unchanged  after  the  anneals.  This  treat¬ 
ment  produced,  in  each  case,  a  p-type  layer  with  a  sheet  resistivity  of  about 
105  Q  /  ^2  »  due  to  a  dominant  acceptor  level  about  0.1  eV  from  the  valence 
band.  Also,  a  1.41-eV  emission  line  was  observed.  However,  we  will  show  that 
the  acceptor  is  definitely  not  related  to  Mn,  even  though  the  1.41-eV  emission 
probably  is.  Thus,  the  elimination  of  Mn  from  the  substrate  and  growth 
ambient  will  not  eliminate  conversion,  at  least  not  in  all  cases.  Indeed, 
it  appears  that  a  defectlike  state  may  be  the  most  important  constituent  of 
the  0.1-eV  center. 
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2.5  The  Effects  of  SI  GaAs  Substrates  on  FET  Devices 
The  following  abstract  is  reprinted  from  Ref.  30. 

A  detailed  investigation  of  21  SI  GaAs  substrates,  Cr-doped,  oxygen- 
doped,  and  undoped,  Bridgman  and  Czochralski,  from  8  different  manufacturers, 
has  been  carried  out  by  temperature-dependent  Hall  effect  (TDH) ,  differential 
Hall  effect  (DH) ,  photoluminescence,  spark-source  mass  spectroscopic  and 
secondary-ion  mass  spectroscopic  measurements.  Also,  Se-ion- implanted 
MESFET's,  with  4  ym  gates,  were  fabricated  on  several  of  the  substrates  and 
tested.  Electrical  conversion  tests  were  carried  out  to  simulate  typical 
VPE  growth  conditions  (750°C,  15  min.,  hydrogen  atmosphere,  no  cap),  and 
ion- implantation  annealing  conditions  (900°C,  15  min.,  hydrogen  atmosphere, 
Si^N^  cap).  The  following  conclusions  have  emerged  from  this  study:  (a)  It 

is  very  difficult  to  get  reliable  n-type  activation  at  less  than  the  mid-10^ 

-3  3  2 

cm  level,  at  least  with  120  Ke  V  Se  ions.  (b)  Low  mobilities  ($10  cm 

V  ^s  ^),  often  seen  in  SI  GaAs  substrate  materials,  are  due  to  inhomogeneity 

g 

(if  p  •<  4  x  10  ft  cm).  The  effective  channel  mobilities  in  implanted 
MESFET's,  near  pinch-off,  are  heavily  influenced  by  these  low  substrate 
mobilities.  (c)  The  commonly  observed  0.5  eV  acceptor  conversion  level  is 
due  to  Fe,  but  the  0.1  eV  acceptor  level  is  not  always  due  to  Mn,  as 
evidently  believed  by  many  workers. 
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3.0 


PHOTOLUMINESCENCE  IN  SEMI-INSULATING  GALLIUM  ARSENIDE 


3.1  Introduction 

A  major  impediment  to  the  development  of  GaAs  integrated  circuits 
and  microwave  devices  has  been  the  lack  of  a  reliable  supply  of  semi-insu¬ 
lating  substrates.  The  main  question  surrounding  semi- insulating  GaAs  has 
concerned  the  compensation  mechanism  by  which  the  undoped  material  is 
semi-insulating.  The  understanding  of  the  compensation  mechanism  has  two 
important  practical  consequences.  First,  knowledge  of  cause-effect 
relationships  between  crystal  growth  and  material  properties  would  greatly 
improve  the  yield  of  semi- insulating  crystals.  Second,  this  understanding 
would  greatly  improve  the  device  performance. 

It  is  well  known  that  the  main  deep  donor  located  at  ~  0.75  eV 
from  the  conduction  band  is  present  in  GaAs  materials  and  that  this  donor 
is  primarily  responsible  for  semi- insulating  properties  through  certain 
compensation  mechanisms.  In  our  study,  deep  center  luminescence  has  been 
studied  through  the  change  of  temperature,  the  persistent  photoluminescence 
quenching  effect,  and  the  stoichiometry  of  crystals,  in  order  to  identify 
the  radiative  mechanisms  involved  in  deep  centers  in  undoped  and  Cr-doped 
semi-insulating  materials. 

3.2  Characteristics  of  Undoped  Semi-insulating  GaAs  (Ref.  31) 

A  photoluminescence  emission  study  has  been  made  on  undoped  LEC 

GaAs  crystals  grown  from  gallium-  and  arsenic-rich  melts.  A  77-meV  deep 
acceptor  responsible  for  the  1.44-eV  emission  is  present  both  in  p-type 
conducting  and  n-type  semi-insulating  crystals.  The  0.68-  and  0.77-eV 
emissions  present  in  the  semi-insulating  materials  are  explained  in  terms 
of  the  radiative  transitions  between  the  band  edge  and  deep  centers. 

Figure  3-1  shows  the  typical  deep-center  photoluminescence  emission  from 
n-type  semi-insulating  GaAs. 
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3.3  Temperature  Dependence  of  0.68-eV  Band  in  GaAs  (Ref.  32) 


The  0.68  eV  photoluminescence  band  present  in  undoped  semi-insulating 
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GaAs  crystals  has  been  studied  with  the  change  of  temperature.  It  is 
shown  that  the  0.68  eV  band  is  due  to  the  radiative  transition  involving 
a  main  deep  donor  and  the  valence  band.  The  origin  of  the  donor  is 
intrinsic  and  may  involve  an  As  antisite  defect.  It  is  found  that  the 
donor  level  does  not  change  in  energy  with  respect  to  the  valence  band  at 
T  =  4-300  K.  The  donor  level  is  found  to  be  at  0.73  eV  from  the  conduction 
band  at  T  =  4K.  The  temperature  dependence  of  the  0.68-eV  band  at  different 
temperatures  is  shown  in  Figure  3-2.  The  temperature  dependence  of  the 
half-width  for  two  samples  is  also  illustrated  in  Figure  3-3. 

3.4  Quenching  Behavior  of  0.68-eV  Band  in  GaAs  (Ref.  33) 

The  persistent  photoluminescence  quenching  effect  is  shown  to  occur 
on  the  broad  0.68-eV  emission  commonly  present  in  undoped  semi-insulating 
GaAs  bulk  materials.  The  dependence  of  this  persistent  emission  on  below 
band-gap  energy  excitation  has  been  measured.  The  data  indicate  involvement 
of  the  main  deep  donor  EL2  in  the  radiative  mechanism  of  the  0.68-eV 
emission.  The  persistent  photoluminescence  quenching  effect  can  be  explained 
by  the  presence  of  both  normal  and  metastable  states  of  EL2.  Figure  3-4 
shows  the  0.68-eV  band  photoluminescence  quenching  rate  vs.  below-bandgap- 
excitation  energy. 

3.5  0.63-eV  Band  Related  to  Oxygen  in  GaAs  (Ref.  34  and  35) 

Temperature-dependent  photoluminescence  and  photoluminescence 

excitation  spectroscopy  have  been  used  to  measure  the  0.63-eV  luminescence 
band  present  in  0-doped  semi-insulating  GaAs.  It  is  shown  that  the  0.63-eV 
band  is  related  to  the  presence  of  0.  The  center  responsible  for  the  band 
forms  a  deep  level  similar  to  the  main  deep  donor  EL2.  However,  spark 
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Fig.  3-3  Temperature  dependence  of  the  half-width  of  the  0.68-eV 
band  for  two  samples. 
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Fig.  3-4 

0.68-eV  band  PL  quenching  rate  vs.  below-band-gap 
excitation  energy.  Solid  line  is  the  representation  of 


w. 


source  mass  spectrometry  indicates  that  incorporation  of  0  into  GaAs  is 
difficult.  A  configuration  diagram  showing  EL2  and  O-related  level  is 
presented  in  Figure  3-5. 

3.6  0.58-eV  Emissions  in  Cr-doped  Semi-Insulating  GaAs  (Ref. 36) 

Heat-treatment  studies  have  been  made  on  n-type  conducting  and 
semi-insulating  GaAs  crystals.  Sharp  zero-phonon  lines  at  0.598,  0.591, 
0.584  and  0.576  eV  and  their  phonon  side-bands  have  been  observed  after 
heat  treatment.  It  is  found  that  the  photoluminescence  emissions  are  due 
to  the  presence  of  Cr.  Close  neighbour  pairs  involving  Cr  and  other 
impurities  are  probably  responsible  for  the  emissions.  Figure  3-6  shows 
the  appearance  of  the  0.58-eV  emission  after  heat  treating  at  600°C.  The 
0.584-eV  zero-phonon  line  and  its  phonon  side-band  are  seen  clearly. 


Lattice  Coordinate 


Fig.  3-5  Configuration  diagram  showing  EL2  and  O-related  level. 
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4  X)  PHOTOLUMI NESCENCE  IN  CONDUCTING  GALLIUM  ARSENIDE  AND  INDIUM  GALLIUM 
ARSENIDE 

4.1  Introduction 

The  incorporation  of  an  intrinsic  defect  has  been  known  to  produce 

electrically  active  centers  in  GaAs.  In  particular,  many  studies  have  been 
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made  on  antisite  defects  formed  during  crystal  growth  *  and  electron  or 
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neutron  irradiation  *  .  The  anion  antisite  double  donor  As„  or  its 

Ga 

complex  has  been  explained  to  be  responsible  for  a  main  electron  trap 
located  at  -0.75  eV  from  the  conduction  band  (commonly  known  at  EL2).  The 
EL2  center  is  a  main  donor  controlling  electrical  properties  in  semi-insulat¬ 
ing  substrate  materials  with  the  compensating  effective-mass  type  shallow 
acceptor  C  ,  as  discussed  in  Section  3.  The  concentration  of  EL2  increases 

AS 

with  the  As-melt  composition.  By  using  the  same  logic  it  can  be  expected 
that  an  acceptor  can  be  produced  by  increasing  the  Ga-melt  composition.  This 
center  is  the  cation  antisite  acceptor  Ga^.  We  have  shown  that  such  a 
center  is  present  in  p-type  conducting  and  n-type  semi- insula ting  GaAs.  We 
will  discuss  the  GaAg  center  and  some  photoluminescence  properties  of  epi¬ 
taxial  InGaAs  crystals. 

4.2  Temperature  Dependence  of  1.44-eV  Emission  in  p-type  GaAs 
(Ref.  40) 

A  study  of  the  photoluminescence  emission  band  at  -1.44  eV  present 
in  GaAs  has  been  made  at  temperatures  between  2-300  K.  Changes  in  photo¬ 
luminescence  excitation  intensity,  emission  energy,  and  emission  intensity 
as  a  function  of  temperature  lead  to  the  identification  of  a  77  +  2-meV 
deep  acceptor.  The  temperature  dependence  of  the  emission-peak  energy 
follows  Eagles'  model  for  the  free  electron-neutral  acceptor  transition. 
Figure  4-1  shows  the  temperature  dependence  of  the  1.44-eV  band. 


Emission  Intensity (Arb.Units) 


Go  As 

LEC  Grown 
p-Type 


-'.  l-  ■V.'S-.T.' 


~l.44  eV 


T=  5  K 


T=  35  K 


T  =  7I  K 


T  =  I09K 


T=I54  K 


1.37  1.39  1.41  1,43  1.45  1.47  1.49  1.51  1.53 

Fig.  4-1  Temperature  dependence  of  the  |T  ^  )  "CAs"  and  "EIn"  resPectively  raean  the  CAg- 

1.44-eV  band  and  the  ^ther  emissions.  26  related  and  intrinsic-region  emissions. 


4.3  Energies  of  the  Double  Acceptor  Ga  in  GaAs  (Ref.  41) 

AS 

Acceptors  present  in  undoped  p-type  conducting  GaAs  have  been 
studied  with  photoluminescence,  temperature-dependent  Hall  measurements, 
deep  level  transient  spectroscopy,  and  spark  source  mass  spectrometry. 

It  is  shown  that  p-type  conduction  is  due  to  the  presence  of  the  shallow 
acceptor  CAg  and  the  cation  antisite  double  acceptor  Ga^.  The  first  and 
second  ionization  energies  determined  for  Ga  are  77  and  230  meV  from 

AS 
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the  valence-band  edge.  Figure  contains  the  4.2-K  photoluminescence 
bands  showing  the  first  and  second  ionization  levels  of  the  Ga^. 

4.4  Photoluminescence  Excitation  of  the  1.44-eV  Emission  in 
p-type  GaAs  (Ref.  42) 

Photolurainescence  excitation  spectra  have  been  obtained  at  4.2  K 
for  the  1.441-eV  cation  antisite  double  acceptor  Ga  emission  present  in 

AS 

undoped  p-type  GaAs.  The  oscillating  structures  at  above-band-gap  energy 
can  be  attributed  to  the  energy  relaxation  of  the  conduction-band  electron 
by  the  successive  emissions  of  the  zone-center  longitudinal  phonon.  The 
energy  of  the  oscillatory  periods  shows  that  the  1.441-eV  emission  at 
4.2  K  is  primarily  due  to  the  neutral  donor-acceptor  pair  transition 
involving  the  Ga  center.  Figure  4-3  shows  the  photoluminescence  excita- 

AS 

tion  spectra  obtained  from  two  samples. 

4.5  Photoluminescence  of  Mn-  and  Un-doped  Ga^  ^In^  53^s  on 
(Ref.  43) 

Photoluminescence  measurements  on  Mn-  and  Un-doped  liquid  phase 
epitaxy  Ga^  ^In^  53^  on  have  been  made  at  the  temperature  range  of 
4.2-296K.  The  Mn  acceptor  binding  energy  is  estimated  to  be  45  +  2  meV. 
Luminescence  due  to  the  Mn  acceptor  quenches  at  60-100K  depending  on  the 
doping  concentration  of  Mn  in  favor  of  a  near-band  edge  emission.  Figure 
4"4  shows  the  4.2-K  photoluminescence  in  undoped  and  Mn-doped  Ing  ^Ga-g  ^ 


GaAs' 

p-type 


Photoluminescence  excitation  spectra  obtained  from  two 
p-type  samples. 


PL  Intensity  (a.u.) 


E(eV) 


Fig.  4-4  4.2-K  photoluminescence  spectra  obtained  from  (a)  undoped 

n-type  samples,  (b)  lightly  Mn-doped  p-type,  and  (c) 
heavily  Mn-doped  p-type  samples. 
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5.0  LOCAL  VIBRATIONAL  MODE  IMPURITY  STUDIES  IN  GALLIUM  ARSENIDE 


Impurity  analysis  of  semiconductor  materials  such  as  GaAs  may  be 
accomplished  with  quantitative  accuracy  and  sensitivity  by  use  of  infrared 
techniques.  Impurities  whose  atomic  masses  are  smaller  than  those  of  the 
host  semiconductor  induce  a  localized  vibrational  mode  (LVM)  absorption 
which  may  be  detected  in  infrared  transmission.  One  important  aspect  of 
such  measurements  is  that  the  location  in  the  lattice  (i.e.,  being  substi¬ 
tutional  on  either  an  As  or  Ga  site)  will  cause  frequency  shifts  of  the 
LVM.  In  addition,  this  technique  is  nondestructive  and  does  not  depend 
on  any  particular  electrical  activation  of  the  impurity  for  it  to  be 
measured. 

The  main  thrust  of  the  investigation  has  involved  high  resolution 
measurements  of  the  LVM  due  to  carbon  and  silicon  in  GaAs.  A  rapid  scan 
Fourier  Transform  Infrared  spectrometer  was  operated  at  the  highest 
attainable  resolution  of  0.06  cm  1  with  the  sample  held  at  80  K  in  a 
cryostat  equipped  with  CdTe  windows.  A  helium  cooled  bolometer  (nitrogen 
cooled  HgCdTe  detector)  with  mylar  (germanium)  beamsplitter  was  used  for 
the  silicon  (carbon)  LVM  measurements.  Samples  with  carbon  concentrations 

ico  -J7  o  18*} 

[C]  between  A  x  10  /cm  to  1  x  10  /cm  and  [Si]  between  1  x  10  /cm  to 
19  3 

5  x  10  /cm  were  found  to  have  additional  structure  as  compared  to  previous 

low  resolution  measurements.  On  the  basis  of  these  measurements,  direct 

evidence  for  the  site  of  substitutional  carbon  in  GaAs  as  being  on  the  As 

44 

sublattice  was  obtained. 

The  cause  of  the  splittings  of  the  LVM  absorptions  for  impurities 
on  the  As  sublattice  are  the  various  configurations  of  the  isotopes  of 
nearest  neighbor,  nn,  Ga  isotopes.  Using  the  relative  natural  abundance 
of  the  Ga  isotopes  of  masses  69  and  71  which  are  60.4%  and  39.6%,  respec¬ 
tively,  and  the  statistics  of  the  possible  configurations  for  the  four  nn 
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G a  Isotopes,  it  was  possible  to  use  a  linear  chain  model  (neglecting  bond 

bending)  to  predict  the  relative  heights  of  the  individual  lines.  Due  to 

several  near  degeneracies,  the  predicted  heights  (assuming  equal  oscillator 

strengths  for  each)  were  1.00,  0.35,  0.24,  and  0.50  whereas  measured 

strengths  in  carbon  doped  samples  was  1.00,  0.4,  0.3,  and  0.5  where  in 

69 

each  case  normalization  was  to  the  highest  energy  or  all  Ga  isotope  nn 
case.  Furthermore,  the  frequency  shift  when  all  nn  are  of  *^Ga  t0  that  of 
all  ^Ga  may  be  calculated  by 


U>. 


/  M  ,  2M„  -  M 
1  \  2  c , Si 


0), 


M  l  2M.  -  M 
2  1  c,Si 
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which  yields  1.0013  for  C  and  1.0036  for  '  Si  as  compared  to  the  measured 


value  of  1.00083  and  1.0028  for  these  impurities.  A  typical  spectrum  for 


carbon-doped  GaAs  is  shown  in  Figure  5-1.  This  particular  sample  was 


17  3 

analyzed  by  Spark  Source  Mass  Spectroscopy  to  have  [C]  =  1.3  x  10  /cm 


17  3 

and  by  Secondary  Ion  Mass  Spectroscopy  to  have  [C]  =  1  x  10  /cm  which 


45  -2 

yields  a  calibration  factor  of  1  cm  absorption  corresponding  to  a  [C] 


of  3.6  x  lO^/cm^. 


One  of  the  most  important  features  of  this  entire  study  is  the 
unambiguous  determination  of  the  site  for  carbon  in  GaAs.  Since  carbon 
and  silicon  can  occupy  either  Ga  or  As  sites,  with  consequent  donor  or 


acceptor  electrical  characteristics,  this  is  of  some  technological  inter- 


46 


est.  In  previous  studies  the  LVM  absorption  due  to  Si  was  identified 


to  occur  at  384  cm  ^  and  that  due  to  Si.  at  399  cm  Such  lines  were 

As 


measured  in  the  high  resolution  condition  to  yield  the  results  shown  in 


Figure  5-2.  The  lineshape  of  the  Si^g  is  extremely  similar  to  that  of 


-1 


carbon  at  582  cm  which  suggests  that  carbon  occupies  the  As  site  and 
that  the  theoretical  explanation  is  correct.  No  line  could  be  identified 


32 


transmission 


»h  resolution  infrared  transmission  spectrum  of 
Licon-doped  GaAs. 


in  carbon-doped  samples  which  would  correspond  to  C  although  the  C 

G3.  As 

13 

line  was  observed.  Since  C  is  only  1%  abundant,  the  CL,  must  be  less 

Ga 

than  this,  if  it  exists  at  all.  This  is  to  be  contrasted  with  the  behavior 

of  Si  in  GaAs,  which  prefers  the  Ga  site  by  two  orders  of  magnitude  over 

the  As  site.  Several  theoretical  calculations  have  now  been  performed 

which  support  the  conclusions  drawn  by  the  simple  linear  chain  model. 

47  48 

Second  order  degenerate  perturbation  theory  and  Green's  function 
techniques  also  fit  the  data  well  and  yield  important  lattice  constant 
information. 

Additional  features  in  Si  heavily-doped  GaAs  were  also  studied 

in  some  detail.  For  example.  Si  -Si.  pair  modes  should  also  show  the 

effects  of  nn  Ga  isotopes.  Since  in  this  case,  the  symmetry  is  lower 

(C^)  and  there  are  only  three  Ga  neighbors  to  the  SiAg  in  the  Si-Si 

pair,  the  lineshapes  are  expected  to  differ  from  the  previous  cases.  The 

high  frequency  mode  (A^+) ,  being  singly  degenerate  splits  into  four  lines 

whose  widths  are  unfortunately  less  than  the  spacing  and  hence  result  in 

an  asymetric  lineshape.  The  low  frequency  mode  at  393  cm  ^  is  doubly 

degenerate  (E+)  which  causes  some  split  lines  to  become  degenerate.  Again, 

the  result  predicted  by  a  linear  chain  model  can  be  made  to  fit  the 

observed  asymetric  lineshape  assuming  only  that  individual  modes  are  of 

typical  widths  as  other  lines  in  the  sample  (0.7  cm  b  with  the  only 

adjustable  parameter  being  the  spacing  between  lines.  Although  the  fit 

is  not  unique,  a  spacing  of  0.4  cm  ^  appears  correct.  A  Green's  function 

49 

calculation  has  been  done  for  the  Si-Si  pairs  which  shows  essentially 


the  observed  behavior.  Finally,  in  a  sample  doped  both  with  Si  and 

28  28  30 

Si,  the  modes  were  resolved  corresponding  to  Si„  -  Si.  and 


^Si  -^Si  to  be  at  455.80  and  457.01  cm  \  respectively,  which  is  again 

03  AS 

a  very  sensitive  measure  of  lattice  spring  constants. 

In  Si-doped  GaAs  which  was  compensated  by  either  ^Li  or  ^Li  dif¬ 


fusion,  the  isolated  Li  line  was  seen  for  the  first  time,  although  it 
had  been  sought  previously  by  other  groups.  This  mode  occurs  at  449.64  cm 
for  ^Li  ,  which  is  near  the  energy  which  might  be  predicted  from  the  \i 


-l 


absorption  in  a  Si  -Li  pair,  when  the  formula 
03  03 


2  .1  5 


3  v  vT 
i-1  1 


is  used  with  the  observed  position  of  454.19,  447.39,  and  437.75  cm  ^  or  a 
predicted  value  of  446.49  cm  A  similar  result  for  ^Li  is  obtained  with 
the  observed  ^Li,,  occurring  at  482.31  cm  Examples  of  such  samples  may 

03 

be  found  in  Figure  5-3. 

In  Si-doped  GaAs  which  was  compensated  with  electron  irradiation, 
additional  structure  was  observed  between  365-376  cm  One  line  in  par¬ 
ticular,  366.63  cm  \  has  been  attributed  to  Si  paired  with  an  nn  defect 

AS 

on  the  Ga  lattice  such  as  the  Ga  vacancy  or  As  antisite"^.  The  large 

03 

width  of  the  line  suggests  Si  involvement.  Also,  the  defect  appears  to 

AS 

have  a  donor  activation  and  appears  in  both  electron-irradiated  and  Li- 
compensated  samples.  Previous  studies^  indicate  that  the  line  diminishes 
with  increased  electron  irradiation  and  grows  upon  annealing.  The  most 
probable  candidate  fitting  these  observations  is  Si  -As  and  is  line 
"A"  in  Figure  5-3.  The  origin  of  lines  "B"  -  "E"  is  undetermined  although 
closely  associated  with  the  electron  irradiations.  An  additional  important 
conclusion  when  comparing  samples  compensated  by  either  method  is  that  in 
high  resolution  neither  the  position  nor  widths  of  any  line  were  different 
between  the  two  compensation  methods.  It  is  well  known,  however,  that 
amplitude  will  vary  due  to  impurity  redistribution.  This  important  result 
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ENERGY  (cm*1) 

Fig.  5-3  High  resolution  infrared  absorption  spectra  for 
two  silicon-doped  GaAs  samples.  Sample  1  was 
compensated  by  electron  irradiation  while 
Sample  2  was  compensated  by  Li  diffusion. 


verifies  that  samples  may  be  studied  using  either  means  of  compensation 
without  perturbation  due  to  the  presence  of  Li.  Additional  unidentified 
absorptions  were  also  investigated  but  no  conclusive  evidence  for  an 
identification  was  obtained. 


6.0  I-V  CHARACTERISTICS  OK  GaAs  JI'ET's 

6.1  Introduction 

Since  the  junction  field-effect  transistor  (JKET)  is  one  of  the 
most  important  devices  made  of  semiconductor  materials  and  the  current- 
voltage  (I-V)  characteristics  are  the  basic  tool  to  study  the  performance 
of  such  devices,  it  is  very  important  to  study  and  understand  the  quanti¬ 
tative  dependence  of  the  I-V  properties  on  the  basic  material  parameters 
obtained  mainly  from  electrical  characterization.  Carrier  drift  velocity 
as  a  function  of  the  electric  field  (V-E)  in  the  conducting  channel  of 
the  device  is  introduced  as  an  input  in  the  calculation  of  the  I-V 
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characteristics.  Previous  models  use  a  V-E  curve  which  is  either  linear 

or  gradually  changes  from  the  low-field  range  of  constant  mobility  to  a 
53 

saturation  value  .  The  latter  V-E  relationship  represents  correctly  the 
behavior  in  semiconductors  like  Si  or  Ge.  However,  for  some  semiconductors 
with  two  conduction-band  valleys,  e.g.,  GaAs  and  InP,  the  velocity-versus- 
field  curve  has  a  peak,  and  a  so-called  negative  conductance  region,  in 
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which  the  drift  velocity  decreases  with  increasing  field.  We  have  shown 
that  a  large  difference  occurs  in  the  I-V  characteristics  and  saturation 
current,  especially  for  small  gate  voltages,  if  one  takes  into  account  the 
realistic  field-dependent  velocity  for  GaAs.  In  the  next  section  we 
briefly  describe  the  model.  A  detailed  discussion  may  be  found  in  Ref.  54. 

6.2  Theoretical  Model  for  JFET 

Lehovec  and  Zuleeg  used  the  Trof imenkof f-type  velocity-field 
relationship  in  their  model 


V  E 

1  +  uE/v 


(6-1) 


Kroemer  introduced  the  following  velocity-field  characteristic 


„r  1  +  (AE) 

v  =  uEt- - r  nr~ 


(6-2) 


1  +  (AE)  pE/v 


m 


At  high  fields,  v  approaches  a  saturation  velocity,  v  and  at  low  fields, 

m 


the  velocity  is  linear  in  E  with  a  low-field  mobility  p.  A  peak  occurs 


in  the  intermediate  region.  This  expression  gives  a  good  fit  to  the 
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measured  velocity  as  shown  in  Figure  6-1,  where  we  have  assumed  the 


saturation  velocity  to  be  1.14  x  10  cm/sec  and  the  low-field  mobility  to 


2  ^ 

be  7500  cm  -V/sec  .  The  values  of  the  two  fitting  parameters  are: 


-1 


A  =  5.2  kV/cm  and  k  =  4.2.  In  the  same  figure,  we  have  also  shown  the 
Trofimenkoff  v-E  curve  (Eq.  6-1)  with  the  same  values  of  low-field  mobil¬ 


ity  p  and  the  saturation  velocity  v  .  We  then  obtain  the  I-V  character- 

m 


istics  using  the  field-dependent  velocity  given  by  Eq.  6-2  and  with  the 
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other  approximations  similar  to  those  used  by  Shockley 

Current  is  proportional  to  the  drift  velocity  v  and  the  open 
channel  thickness  (a-d) ,  where  a  and  d  are  the  total  channel  and  depletion 
thicknesses,  respectively.  Since  the  Poisson's  equation  relates  this  d 
to  the  potential  at  that  point  and  v  is  a  function  of  the  field,  which 
is  just  the  derivative  of  the  potential,  we  obtain  a  differential  equation 
for  the  potential  V,  with  the  current  as  a  parameter  in  the  equation. 
Solving  the  differential  equation  into  proper  boundary  conditions,  namely. 


X  -  0,  ,  -  t  =  <-VCB/Vo) ! 


(6-3) 


x  -  L,  y  -  u  E  [<VD-VCB/VoP  , 


where  V  is  the  sum  of  the  gate  voltage  and  the  built-in  voltage,  V  ,  is 
G15  D 


) 


the  applied  drain  voltage,  x  =  0  and  L  are  the  positions  of  the  source  and 

drain,  respectively,  and  is  a  parameter  (which  depends  on  impurity 

density  in  the  channel  and  the  total  channel  thickness)  we  obtain  an  implicit 

equation  for  the  current  as  a  function  of  V  ,  V  and  their  material  and 

JD  O 

device  parameters.  With  the  V-E  relationship  given  by  Eq.  6-2,  it  is  not 

possible  to  solve  this  equation  analytically,  although  its  numerical 

solution  is  simple.  Figure  6-2  shows  the  numerical  solution  for  the 

reduced  current  3  with  the  reduced  gate  voltage  t  as  a  parameter  (solid 

curves).  The  dashed  curves  are  plots  of  the  results  of  Lehovec  and  Zuleeg 
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corresponding  to  the  velocity  characteristics  represented  by  the  dashed 

line  in  Figure  6-1.  In  all  of  the  figures  solid  lines  are  referred  to  as 

CL  and  dashed  ones  as  LZ.  The  various  parameters  used  in  our  calculations 

are:  N  =  3  x  1015cm  W  =  1.25  mm*  a  =  1.4  pm,  L  =  25  ym,  e  =  12.9.  These 

53 

are  the  parameters  for  the  Device  PN-l-S  of  Lehovec  and  Zuleeg  .  Conclu¬ 
sions  drawn  in  this  paper  also  hold  for  devices  with  more  realistic  para¬ 
meters  for  GaAs  FET's.  For  a  better  understanding  of  the  reasons  for  the 
large  difference  in  the  I-V  curves  (solid  and  dashed)  we  have  computed  the 
open  channel  thickness  (i.e.,  the  shape)  and  the  field  along  the  channel  and 
have  plotted  these  quantities  in  Figures  6-3  and  6-4,  respectively.  We  have 
plotted  in  Figure  6-5  the  transconductance  (in  units  of  the  transconduc¬ 
tance  g^o  for  the  Shockley  model  (i.e.,  linear  v-E  relationship)  versus 
the  reduced  gate  voltage ' -V„D/V  .  The  transconductance  obtained  by  Lehovec 

CjD  O 
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and  Zuleeg  does  not  agree  with  that  of  the  present  calculation.  This 
shows  that  the  CL  saturation  currents  for  various  gate  voltages  cannot  be 
obtained  by  a  simple  scaling  from  those  of  LZ. 


Fractional  thickness  of  the  open  channel,  as  a 
function  of  position,  for  the  parameters  shown 


6.3  Physical  Parameters  of  FET's 

In  the  model  described  in  Section  6-2,  we  have  fitted  the  experi¬ 
mental  V-E  data  of  Smith,  et  al.^  with  Eq.  6-2.  This  set  of  data  is  for  a 

2 

pure  sample,  whose  low-field  mobility  (p  )  is  7500  cm  /V-sec.  However, 

Li 

the  conducting  channel  region  of  an  actual  device  has  lower  p  .  The  main 

L 

reason  for  this  low  p  is  that  impurities  are  incorporated  in  the  region 

Li 

to  provide  the  carriers  for  the  device  operation.  Scattering  of  the 

carriers  from  the  ionized  impurities  reduces  the  mobility  at  lower  fields. 

Since  the  Coulomb  scattering  cross  section  is  inversely  proportional  to 

the  fourth  power  of  the  velocity  of  the  incident  particle,  the  V-E  curve 

is  not  expected  to  be  much  different  at  higher  fields.  Therefore,  in  the 

absence  of  measured  V-E  curve  for  the  channel  material,  we  need  to  modify 

the  V-E  relationshp  (Eq.  6-2)  such  that  it  gives  p  equal  to  the  measured 

Li 

2 

value  (-.4000  cm  /V-sec),  but  it  remains  the  same  in  the  high  field  region. 

2 

To  achieve  this,  we  add  an  emperical  term  (<*E  )  such  that  the  peak  in  V-E 
remains  the  same.  We  have  developed  an  efficient  computer  program  to  fit 
the  measured  I-V  data  with  the  calculated  I-V  data,  and  the  best  fit  pro¬ 
vides  the  geometrical  and  material  parameters.  The  model  calculation 
includes  the  aforementioned  modified  V-E  relationship.  It  may  be  noted 
that  if  a  more  reliable  V-E  relationship,  either  measured  or  theoretically 
calculated,  is  available,  we  would  be  able  to  incorporate  it  fairly  simply. 

The  following  are  the  parameter  values  for  the  device  COM  130S87 
obtained  by  fitting  the  I-V  data  (measured  in  this  laboratory)  with  our 
model  (See  Figure  6-6). 

total  channel  thickness  =  0.25  p 

gate  length  =  9.8  p 


built-in  voltage 


-1.0  v 


With  these  parameter  values  we  obtain  a  very  good  fit.  The  LZ  model  gives 
a  fit  five  times  worse  than  our  model.  Considering  the  uncertainties  in  the 
V-E  curve  and  that  there  is  no  unphysical  parameter  in  the  model,  we  may  say 
that  the  parameters  obtained  and  the  fit  are  quite  good. 

The  main  reason  that  the  fit  gets  worse  at  high  gate  and  drain 

voltages  (V  and  V  )  is  that  for  high  V  and  V  there  is  some  conduction 
li  U  Cj  D 

through  the  substrate  and  consequently  the  current  (1^),  instead  of  reaching 
a  saturation,  has  a  small  slope.  This  factor  is  not  incorporated  in  the 
model. 

6.4  Conclusion 

In  conclusion,  we  have  shown  that  for  FET's  composed  of  semicon¬ 
ductors  like  GaAs  and  InP,  the  theoretical  I-V  characteristics  can  be  in 
error  by  a  large  amount  if  one  uses  the  LZ  model,  even  though  the  velocity 
(Eq.  6-1)  is  correct  at  very  small  and  very  high  fields.  This  error  arises 
because  the  Trofimenkoff  relationship  differs  from  the  realistic  v-E  char¬ 
acteristic  for  the  materials  under  consideration  in  the  intermediate  field 
region.  Using  a  realistic  v-E  relationship  we  have  developed  a  model  which 
gives  a  good  fit  with  measured  I-V  data.  We  hope  that  this  model  can  be 
modified  and  used  in  analyzing  the  I-V  data  of  HEMT  devices  which  are  under 
study  in  the  current  program. 


7.0  QUANTUM  WELI-S  AND  SUlMiKLATTICES 


7.1  Introduction 

Ln  recent  years  an  intense  worldwide  effort  has  been  devoted  to 
the  studies  of  the  properties  of  a  new  class  of  man-made  materials,  called 
superlattices.  Semiconductor  superlattices  are  formed  by  growing  alternate 
layers  of  two  different  lattice  matched  (or  slightly  mismatched)  semicon¬ 
ductors.  The  most  studied  system  is  that  of  GaAs  and  Ga,  A Z  As.  Growth 
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techniques  such  as  Molecular  Beam  Epitaxy  (MBE)  or  Metal  Organic  Chemical 
Vapor  Deposition  (MOCVD)  make  it  possible  to  precisely  control  concentra¬ 
tion  of  the  components  (e.g.,  of  A Z  in  GaA£As)  and  to  produce  sharp 
interfaces.  These  systems  have  important  practical  implications.  Owing 
to  the  difference  in  the  band  gaps  of  GaAs  and  GaAftAs,  carriers  (electrons) 
from  the  intentionally  doped  GaA£As  region  tunnel  to  the  lower  energy  GaAs 
region.  Thus,  the  carriers  are  geometrically  separated  from  the  impurity 
ions  and  thereby  the  transverse  mobility  increases  tremendously.  Exploit¬ 
ing  this  high  mobility  and  spatial  quantization,  a  number  of  novel  devices 
have  been  made,  including  high  electron  mobility  transistors  (HEMT)  and 
lasers . 

If  the  GaAs  layers  are  separated  by  thick  layers  of  GaAJtAs 
(>200  A°) ,  the  carriers  in  one  well  region  (GaAs)  do  not  feel  the  presence 
of  other  GaAs  regions.  Then  each  GaA£As-GaAs  -  GaA£As  period  can  be 
thought  of  as  an  isolated  entity,  called  a  quantum  well.  Due  to  the 
quasi-two-dimensional  nature  of  the  carriers  in  such  quantum  wells,  the 
carriers  exhibit  various  interesting  properties  which  are  of  fundamental 
nature.  Thus,  studies  related  to  quantum  wells  and  superlattices  are 
important  not  only  from  a  practical  point  of  view,  but  also  from  a  point 
of  view  of  basic  physics. 


As  in  bulk  semiconductors,  understanding  the  properties  and  per¬ 
formances  of  devices  made  of  superlattices  depends  crucially  on  the  knowl¬ 
edge  of  energy  levels,  effect  of  external  perturbations  on  these  levels, 
quality  of  the  interfaces,  etc.  In  the  following  sections  we  briefly 
describe  some  of  the  contributions  we  have  made  in  this  area. 

7.2  Hydrogenic  Impurity  Ground  State  in  GaAs-GaA£As  Multiple- 
Quantum-Well  Structures 

The  following  abstract  is  reprinted  from  Ref.  57: 

We  report  the  results  of  a  variational  calculation  for  the  hydro- 
genic-irapurity  ground  state  in  a  multiple-quantum-well  structure  consisting 
of  alternating  slabs  of  GaAs  and  Ga^^A&^As.  Calculations  have  been  carried 
out  with  the  assumption  that  the  impurity  envelope  wave  function  spreading 
beyond  the  next-nearest-neighbor  GaAs  wells  is  negligible.  Impurity 
envelope  wave  functions  have  been  plotted  for  some  typical  GaAs  well  and 
Ga^_xA2.  As  barrier  thicknesses  to  find  the  extent  of  wave-function  spreading 
The  binding  energy  is  found  to  vary  substantially  as  a  function  of  the 
barrier  thickness.  Calculations  are  performed  for  the  variations  of  the 
binding  energy  as  a  function  of  the  well  thickness  and  also  as  a  function  of 
the  barrier  thickness.  The  main  peak  in  the  impurity  binding  energy  in 
superlattices  with  equal  well  and  barrier  thickness  is  shifted  towards  a 
thickness  larger  than  that  in  single-well  systems.  A  secondary  peak  appears 
at  a  very  small  thickness,  which  arises  because  the  model  includes  only 
three  wells.  The  results  of  the  present  calculation  in  various  limiting 
cases  agree  with  previous  results. 

The  introduction  from  the  same  paper  is  also  reprinted. 

The  recent  advances  in  crystal-growth  techniques  such  as  molecular- 
beam  epitaxy  and  metal-organic  chemical-vapor  deposition  (MOCVD)  have  made 
possible  the  growth  of  systems  consisting  of  alternate  layers  of  two 


different  lattice-matched  semiconductors  with  very  precisely  controlled 
thicknesses  and  sharp  interfaces.  Such  one-dimensional  periodic  struc¬ 
tures  are  generally  referred  to  as  superlattices.  Numerous  studies  have 

been  devoted  to  various  aspects  of  the  electronic  states  associated 
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with  such  systems  .  Among  the  superlattices  grown,  so  far,  the 
GaAs-Ga1  A£  As  system  is  the  simplest  and  the  most  extensively  studied. 

I'-X  X 

The  band  gap  of  Ga,  A£  As  increases  with  the  concentration  of  x  of 
aluminum.  Thus  there  exists  a  discontinuity  of  the  band  edges  at  the 
GaAs-Ga^^AJl^As  interfaces.  The  conduction-band-edge  discontinuity  is 
about  85%  of  AE  ;  the  band-gap  difference  between  bulk  Ga,  Ail  As  and 
GaAs.  The  valence-band-edge  discontinuity  is  about  15%  of  AE  .  As  a 

s 

result,  the  electrons  and  holes  in  a  GaAs  layer  find  themselves  in 

approximately  rectangular  potential  wells  for  sharp  interfaces. 

Studies  of  the  shallow  impurity  states  in  such  quantum  wells 
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have  recently  begun  with  the  work  of  Bastard  .  He  calculated  the 

ground-state  binding  energy  of  a  hydrogenic  shallow  impurity  in  the 

GaAs  quantum  well.  In  that  calculation  the  potential  barrier  height 

in  the  Ga^_xA£^As  regions  was  taken  to  be  infinite.  He  finds  that  the 

binding  energy  increases  as  the  well  size  (i.e.,  the  GaAs  layer  thickness) 
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is  reduced.  Mailhiot,  Chang,  and  McGill  (MCM)  have  done  an  extensive 
calculation  for  the  binding  energies,  wave  functions,  and  their  varia¬ 
tions  with  well  thickness,  impurity  position,  etc.,  for  realistic  finite 
potential-barrier  height,  contrary  to  Bastard’s  infinite  barrier  height. 
They  also  include  the  effects  due  to  different  effective  electronic 
masses  and  dielectric  constants  in  GaAs  and  Ga^  xA£xAs  layers.  They 
have  found  that  the  binding  energy  goes  through  a  maximum  as  the  well 
size  is  reduced  instead  of  continuously  increasing  as  is  found  in  the 


infinite-barrier  calculation.  Greene  and  Bajaj  (GB)  have  also  performed 
a  similar  calculation  for  impurities  at  the  center  of  a  quantum  well.  Both 
of  the  groups,  MCM  and  GB,  have  calculated  some  excited-state  binding 
energies. 

In  all  of  the  previous  calculations  it  has  been  assumed  that  the 

Ga^_xA£xAs  layers  are  thick  enough  to  confine  the  wave  functions  so  that 

they  do  not  spill  over  to  the  adjacent  GaAs  quantum  wells.  Calculations 

were  performed  with  this  assumption  for  an  impurity  in  a  quantum  well  with 

two  semi-infinite  Ga,  A£  As  barriers  on  each  side  of  the  well.  Super- 
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lattices  are  made  with  layer  thickness  ranging  from  a  few  monolayers  to 

about  400  X.  Most  attention  has  been  focused  on  systems  with  aluminum 

concentration  x  of  Ga^_xA£xAs  less  than  0.45.  In  this  concentration  range 
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the  band  gap  is  direct  at  the  T  point  .  The  spreading  of  the  impurity 

envelope  wave  functions  depends  on  the  potential  barrier  height  as  well  as 

the  barrier  thickness.  In  general  wave  functions  spread  more  to  the 

adjacent  wells  if  the  barrier  height  or  the  thickness  is  reduced.  Thus 

the  previous  calculations  with  single-well  approximation  are  not  adequate 

for  thin  superlattices,  or  even  for  moderately  thick  superlattices  but  with 

small  aluminum  concentrations.  This  is  evident  from  the  wave-function  plot 
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of  MCM  .  In  this  situation  we  feel  that  it  is  desirable  to  perform  a 
calculation  which  should  be  valid  for  thin  or  small  barrier-height  super¬ 
lattices.  The  present  paper,  to  our  knowledge,  is  the  first  report  of  such 
a  calculation.  In  this  calculation  we  consider  for  simplicity  the  impurity 
to  be  at  the  center  of  a  well.  Assuming  that  the  spread  of  the  wave 
functions  to  the  next-nearest-neighbor  wells  of  the  one  containing  the 
impurity,  the  calculation  has  been  performed  with  only  one  well  on  each  side 


of  the  well  under  consideration 


One  of  the  main  features  of  this  calculation  is  presented  in 
Figure  7-1.  The  reader  is  referred  to  the  original  Reference  57  for 
details. 

7.3  Effect  of  Nonparabolicity  on  the  Energy  Levels  of 
Hydrogenic  Donors  and  Excitons 

The  following  abstract  is  reprinted  from  Reference  62: 

Binding  energies  of  the  ground  state  and  of  a  few  low-lying 
excited  states  of  a  hydrogenic  donor  in  a  quantum-well  structure  consist¬ 
ing  of  a  single  layer  of  GaAs  sandwiched  between  two  semi-infinite  layers 
of  Ga^^AJl^As  are  calculated,  including  the  effect  of  nonparabolicity  of 
the  conduction  band  and  following  a  variational  approach.  The  effect  of 
nonparabolicity  of  the  conduction  band  is  included  by  using  an  expression 
for  the  energy-dependent  effective  mass  based  on  the  k-p  approximation. 
The  variations  of  the  binding  energies  of  these  states  as  a  function  of 
the  size  of  the  GaAs  quantum  well  for  different  values  of  the  potential 
barrier  (or  equivalently  for  different  values  of  A£  concentration  x)  are 
calculated.  These  results  are  compared  with  those  obtained  with  the  use 
of  a  parabolic  conduction  band. 

The  following  abstract  is  reprinted  from  Reference  63: 

Binding  energies  of  the  ground  states  of  Wannier  excitons  in  a 
quantum  well  structure  consisting  of  a  single  layer  of  GaAs  sandwiched 
between  two  semi-infinite  layers  of  Ga^  ^AJl^As  are  calculated  including 
the  effects  of  non-parabolicity  of  the  conduction  band,  following  a 
variational  approach.  Due  to  reduction  in  symmetry  along  the  axis  of 
growth  of  these  structures  and  the  presence  of  energy  band  discontinui¬ 
ties  at  the  interfaces,  the  degeneracy  of  the  valence  band  of  GaAs  is 
removed,  leading  to  the  formation  of  two  types  of  excitons,  namely,  the 
heavy  hole  exciton  and  the  light  hole  exciton.  The  effect  of 
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non-parabolicity  of  the  conduction  band  is  included  by  using  an  expression 

— v-  -f 

for  the  energy  dependent  effective  mass  based  on  k  -  p  approximation.  The 
variations  of  the  binding  energies  of  these  two  excitons  as  a  function  of 
the  size  of  the  CaAs  quantum  wei 1  for  various  values  of  the  heights  of  the 
potential  barrier  are  calculated.  These  results  are  compared  with  those 
obtained  using  a  parabolic  conduction  band.  It  is  found  that  the  inclusion 
of  the  non-parabolicity  effect  leads  to  more  binding  for  all  values  of  the 
quantum  well  size. 

7.4  Effect  of  Interface  Quality  on  Photoluminescence  Lineshape 

The  following  abstract  is  reprinted  from  Reference  64: 

We  have  developed  a  simple  theory  to  understand  the  role  of  inter¬ 
facial  quality  in  the  lineshape  of  photoluminescence  spectra  in  quantum 
wells.  The  interface  is  described  in  terms  of  microscopic  fluctuations 
iS^  and  where  6^  is  the  local  fluctuation  in  the  well  width  and  6  2  is 
the  lateral  correlated  extent  of  the  fluctuation.  We  make  use  of  Lifshitz 
theory  of  disordered  alloys  to  determine  the  probability  of  distribution 
of  fluctuations  in  the  well  size  over  the  extent  of  the  optical  probe, 
i.e.,  the  exciton.  The  lineshape  is  then  calculated  from  this  probability 
distribution.  Both  6^  and  are  found  to  be  important  in  controlling 
the  linewidtbs  in  quantum  wells.  The  use  of  this  quantitative  theory  to 
characterize  the  microscopic  nature  of  interfaces  is  discussed. 


8.0  SPARK  SOURCE  MASS  SPECTROMETRY 

8.1  Background 

Spark  source  mass  spectrometry  has  been  utilized  for  over  two 
decades  for  the  comprehensive  elemental  characterization  of  inorganic 
solid  materials,  particularly  with  regard  to  trace  and  ultra-trace  impuri 
ties.  The  technique  features  absolute  qualitative  identification  of 
elemental  components,  with  rapid,  simultaneous  detection  of  most  elements 
to  0.1  parts  per  million,  atomic,  or  lower.  The  role  of  the  SSMS 
(MS-702,  AEI-CED,  Ltd.,  Manchester,  England)  facility  located  at  Wright 
State  University  has  been  to  assist  AFWAL/AADR  in  the  analysis  and 
evaluation  of  gallium  arsenide  research  materials.  Since  the  type  and 
concentration  of  many  elemental  impurities  -  both  doped  and  residual  -  ar 
known  to  drastically  affect  the  electrical  properties  of  GaAs,  the  basic 
function  of  the  SSMS  effort  has  been  to  provide  routine,  reliable  deter¬ 
minations  of  low-level  impurities  in  GaAs  samples.  During  the  past  three 
year.-.,  approximately  200  GaAs  samples,  including  thick  epitaxial  layers, 
doped  samples,  and  substrate  materials,  have  been  analyzed. 

8.2  Analytical  Program 

The  SSMS  technique  has  been  employed  for  three  primary  applica¬ 
tions  : 

1)  Routine  analysis  of  high-purity  GaAs  for  residual  elemental 
impuri ties . 

2)  The  continuation  of  calibration  studies  to  quantify  measure¬ 
ments  of  elements  of  interest  such  as  Fe,  Mn,  Cr,  S,  Si  and 

8. 

3)  The  development  of  methodology  for  determination  of  trace 
levels  of  oxygen  and  carbon. 


In  pursuing  the  first  goal,  the  spark  source  mass  spectrometr ic 
technique  has  been  shown  to  be  effective  in  the  rapid  determination  of  a 
wide  range  of  elemental  impurities.  As  is  shown  in  Table  8-1,  detectabil¬ 
ities  are  frequently  less  than  <0.1  parts  per  million,  atomic  (ppma). 

During  the  three-year  investigation,  many  CaAs  samples  were  found  to  con¬ 
tain  residual  impurities  -  particularly  Fe  and  Si  -  which  (as  revealed 
by  Hal  1-effect  tests)  adversely  affected  the  resistivity  properties. 

In  conjunction  with  the  second  goal,  CaAs  wafers  were  implanted 
using  3V.  Additionally,  some  CaAs  samples  were  implanted  using 
The  results  obtained  from  SSMS  analyses  are  shown  in  Table  8-2.  This  data 
shows  the  SSMS  method  to  be  effective  in  tracking  implanted  doses  and,  in 
the  case  of  3(>0+  a  calibration  factor  ("True"  value/Observed  value)  of 
2-3. 

The  third  goal  sought  during  the  program  has  involved  the  utiliza¬ 
tion  of  a  cryogenic  pumping  system^3  developed  for  the  MS-702  at  Wright 
State  University.  A  machine  modification  involved  the  installation  of  a 
turhomo 1  ecu  I  a r  vacuum  pump  (Alcatel  Vacuum  Products,  TMP  140)  for  the  ion 
source  chamber.  The  purpose  was  to  achieve  "cleaner"  pumping  (reduced 
hydrocarbon  background)  than  was  provided  by  the  standard  oil  diffusion 
pump,  thereby  allowing  lower  limits  of  detection  for  carbon.  The  addition 
of  turbo  pumping  along  with  the  liquid  helium  cryopumping  of  the  ion  source 
region  has  yielded  detectabilities  of  0.02  ppma  for  oxygen  and  0.02  -  0.03 
ppma  for  carbon.  Along  with  the  ion-implantat ion  experiments  described 
above,  (la As  samples  (and  one  Si  sample)  doped  with  oxygen  or  carbon  have 
been  analyzed  and  these  efforts  are  summarized  in  Table  8-3. 
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TABLE  8-1.  SPARK  SOURCE  MASS  SPECTROSCOPIC  DATA  IN  GaAs 


Units:  Parts  per  million  atomic  (ppma) 
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liquid  helium 


8.3  Summary 


Throughout  the  program,  the  application  of  spark  source  mass 
spectrometry  has  been  shown  to  be  generally  successful  in  surveying  GaAs 
and  other  semiconductors  ror  trace  impurity  content.  However,  the  need 
lor  continued  calibration  efforts  using  suitably  doped  samples  is  clear. 
During  the  program,  several  GaAs,  Ini’,  and  GaP  materials  that  had  been 
heavily  doped  with  elements  such  as  Fe,  Mn,  Cr,  Si,  Zn,  Te,  Sn,  and  Ge 
were  analyzed  by  SSMS  and  the  acquired  values  used  in  conjunction  with 
photoluminescence  measurements.  However,  insufficient  sample  quantity 
prevented  analysis  by  complimentary  analytical  techniques.  The  results 
in  Table  8-3  indicate  reasonably  good  SSMS  accuracy  for  0  and  C  deter¬ 
mination  when  suitable  samples  are  available,  and  liquid  helium  cooling 
and  turbomolecular  pumping  are  used.  Subsequent  SSMS  analyses  of  GaAs 
samples  thought  to  contain  0  (as  shown  by  photoluminescence  spectra)  did 
not  reveal  significant  0  levels.  This  could  be  due  to  greater  sensitivity 
by  photo  l  uminescence ,  in  observing  0  peaks,  than  SSMS. 
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9.0  NEW  PAIRING  STATE,  CHARGE  DENSITY  WAVE,  SPIN  DENSITY  WAVE,  PAIRING 
OK  HOLES,  AND  SUPERCONDUCTIVITY 

Tin*  following  summary  is  from  an  invited  talk  given  by  S.  B.  Nam 
at  the  9tli  International  Conference  on  High  Pressure,  in  1983: 

We  have  suggested  that  the  pairing  of  holes^,  or  finite  momentum 
pairing(q),  or  charge  density  wave  (CDW)  orders  might  be  responsible  for 


the  unusual  metastable  phenomena  observed  in  pressure  quenched  CdS 


67,68,69 


Furthermore,  it  is  shown  that  the  pair  binding  energy  in  one  dimension^ 

and  two  dimensions^,  has  an  extra  maximum  at  the  pair  momentum, 

q0  =  k  +  K  ,  additional  to  the  usual  one  at  q  =  0,  where  K  =  (K  +  2mw  ) 
c.  r  r  r  r  C 

and  u)t<  is  the  usual  cut-off  energy  (h  =  1  =  c).  These  findings  encourage 

us  to  believe  that  there  is  a  possibility  of  having  the  finite  momentum 

pairing  order.  This  is  what  we  have  found. 


The  coexistence  of  CDW  and  superconductivity  has  been  examined  by 
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several  people,  in  one  dimension  ,  and  in  A15  compound  such  as  V^Si  , 
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and  others.  Recently  the  observed  Raman  lines  in  2H-NbSe2>  can  be 
understood  by  the  CDW  induced  modes  via  coupling  between  the  CDW  and  the 
superconducting  gap  excitation,  and  in  a  quite  different  view,  also  by 
a  collective  mode  based  on  the  gauge  invariant  theory^.  The  S0(6)  has 
been  identified  as  a  dynamical  group  of  a  mean  field  Hamiltonian  of  the 
CDW  superconductor^. 

On  the  other  hand,  Machida,  et  al.^  have  suggested  the  finite 

momentum  pairing  in  antiferromagnetic  superconductors.  However,  Nass, 
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et  al.  have  pointed  out  that  the  calculation  of  Machida,  et  al.  was 
in  error,  and  have  argued  that  within  their  numerical  calculation,  finite 


momentum  pairing  is  possible  in  neither  antiferromagnetic'  y  nor  CDW 


superconductor 


But,  both  groups  of  Machida,  et  al.  and  Nass,  et  al.  have  treated 
tlu*  order  parameter  as  real  variables.  We  believe  that,  physically,  the 
order  parameters  should  be  treated  as  complex  variables  when  both  zero  and 
finite  momentum  pairing  orders  are  present,  since  the  latter  is  equivalent 
to  having  a  local  current,  even  though  the  total  current  is  zero,  with  the 
q  and  -q  pairings. 

In  fact,  we  have  found  that  the  lowest  free  energy  state  is  indeed 
the  one  in  which  the  phases  of  order  parameters  have  special  values.  This 
is  the  most  striking  result,  in  our  view.  The  state  with  real  order  par¬ 
ameters  turns  out  not  to  be  the  lowest  free  energy  state,  when  the  q  order 
is  present  with  the  other  orders,  such  as  the  BCS,  or  both  CDW  and  BCS 
orders . 

We  have  shown  that  when  the  q  order  is  present  with  the  other 
orders,  the  phases  of  order  parameters  play  crucial  roles  for  determining 
the  ground  state,  and  that  the  excitation  energy  and  effective  energy  gap 


10.0 


ORIGIN  OK  THE  X-CENTERS  IN  SILICON  AND  GERMANIUM;  SOLITON 


For  the  last  several  years,  there  have  been  reported  the  unknown 
gQ_ 33  84~85 

X-centers  in  doped  silicon  ,  and  doped  germanium  .  Several 

attempts  with  models  of  complexes  have  been  made  to  understand  the  X- 

80  82 

centers  with  qualitative  success  *  .  It  has  been  an  important  problem 

to  understand  the  X-centers  since  the  above  materials  have  been  used  for 

sensitive  IR  detectors.  The  materials  with  the  X-centers  cannot  be  used 

for  such  detectors,  since  an  unwanted  noise  is  produced  via  the  X-centers. 
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Recently,  Ohmer  and  Lang  have  been  able  to  obtain  a  good  agree¬ 
ment  between  the  calculated  and  observed  X-center  energies  for  the  majority 
of  data,  using  an  emperical  relation: 


Ex  -  <  ■  *EA  +  V  -  «EA  +  EB> 


(10-1) 


where  E  ,  E  ,  and  E  are  the  energies  of  the  X-center,  and  of  two  appro- 
x  A  B 

priate  acceptors  A  and  B,  respectively.  The  superscript  G  stands  for  the 
respective  ground  states. 

Wo  present  here  a  physical  picture  leading  to  the  above  emperical 
relation  by  demonstrating  a  possibility  of  having  the  distortion- induced 
soli  ton.  By  no  means  is  our  description  a  complete  answer.  Perhaps,  it 
will  give  a  direction  by  which  the  understanding  of  the  X-centers  will  be 
enriched. 

Our  physical  scenario  goes  as  follows:  When  a  distortion  is 
introduced  in  the  system  via  doping,  or  irradiation,  or  other  means,  the 
probability  of  transferring  a  hole  from  the  acceptors  A  to  B,  or  vice 
versa,  will  be  increased.  Of  course,  the  same  argument  will,  equally  well, 
apply  to  the  case  of  electrons  from  two  appropriate  donors.  If  the  transfer- 
integral  energy  exceeds  the  energy  difference  between  the  levels  of  two 


acceptors,  ',.(K^  -  E  ),  then  a  iK*w  order  ground  stale  sols  in,  and  a  solilon 
state  is  predicted  to  occur  with  the  energy  of  k(K  +  H  ) .  This  yields  the 
empirical  relation  of  Eq.  10-1. 

The  two  appropriate  states,  of  course,  are  the  same  symmetry  states 
from  acceptors  A  and  B,  since  other  combinations  will  have  negligibly  small 
or  zero  overlapping  integrals  of  wave  functions.  We  thus  discuss  the  two 
same  symmetry  states,  from  two  appropriate  acceptors,  for  our  purposes. 

The  Hamiltonian  of  the  system  can  be  written  with  4>+  =  (C+  ,  C~*")  as 

A  B 


H  =  Z 


,  + 

cf>  G 


<(>  +  D  +  H 


hD 


(10-2) 


hi) 


-  2  g: 


°  C1  CB  +  C' 


c. 


where  and  (i  =  A,  B)  are  the  usual  creation  and  annihilation  operators 

for  holes,  and  CK  (i  =  1,2,3)  are  the  2x2  Pauli  matrices.  The  second  term 

is  the  distortion  energy  and  is  the  transfer  coupling  energy  between  the 

distortion  field  I)  and  holes,  with  the  transfer  integral  energy  g.  The  e 

is  'j(E  -  K  ),  and  the  zero  energy  is  chosen  at  !2(E.  +  E_)  ,  with  the  positive 
A  d  A  B 

energy  axis  opposite  to  that  of  the  electron.  The  summation  should  be  under¬ 
stood  in  the  usual  manner  over  the  up  and  down  states  of  holes,  for  simplic¬ 
ity.  The  distortion  field  can  be  from  the  lattice  distortion,  or  the  dis¬ 
tortion  of  acceptor  positions  from  their  equilibrium  positions. 

Using  the  usual  commutation  of  CA  with  H,  one  gets  the  eigenvalue 
equat ion , 

(  Uj  c  +  Oj  A)  4>  =  E  <{>  (10-3) 

i 

with  A  =  <-  g 2  D 


where  <  >  denotes  the  grand  canonical  ensemble  average. 


We  easily  got  the  eigenvalues  +  H  and  the  spinor  solution  (J)  =  (  u,  v  )  as 

E  =  (  c2  +  A2  j'2  (10-4) 

u2  =  1  -  v2  =  (  1  +  £/E)/2 

with  sign  (v)  =  sign  (A).  The  Hamiltonian  Eq.  10-2  can  be  rewritten,  in 

+  +  +  + 

terms  of  new  Fermion  operators  (a  ,  a„)=(C.  ,C_)(u+  io„  v  ),  and 

1  Z  Ad  Z 

ni  ”  at  "j  (i  =  1,2)  as 

H  =  X  E(nL  -  n2)  +  D2  (10-5) 

In  the  ground  state,  |I1  a^  >  ,  we  get  n^  =  0,  and  n^  =  1,  so  that  the  new 
up  and  down  states  are  occupied.  The  ground  state  energy  W(D)  is  given 
as 

W(D)  -  W(0)  =  2(e  -  E)  +  D2  (10-6) 

This  has  a  maximum  at  D  =  0,  and  a  minimum  value 

W(  +  I)  )  -  W(0)  =  -  E(  1  -  e/E)2  (10-7) 

—  m 

at  D  =  +  as  shown  in  Figure  10- l,  with  the  condition 

2  2  2  2 

A  »  «  -  C  -  g  D*.  (10-8) 

The  above  condition  is  essential  for  our  discussion.  In  order  to  have 
an  order  parameter  A,  the  c  should  be  less  than  g.  When  it  is  such  a 
new  order  ground  state  sets  in,  so  that  the  distorted  state  (D/0) 
becomes  the  lowest  energy  state,  in  analogy  to  a  charge  density  wave. 


or  the  theory  of  superconductivity.  Note  that  there  are  double  minima 


10-1 


The*  ground  state  energy  as  a  function  of  D  is  in 
the  upper  graph.  The  order  parameter  A,  soliton 
formation  energy  E^,  and  the  minimum  ground  state 

energy  as  functions  of  energy  difference  C  are 

shown  in  the  lower  graph. 

L.  Q 


I)  =  +  I)  .  II"  wo  consider  the  distortion  of  distani'o  between  two  appropriate 
acceptors  for  Llio  disLorLion  field,  then  there  are  two  separation  distances 
with  the  same  energy  of  the  system.  This  double  degeneracy  leads  to  a 
soliton  state.  In  other  words,  one  configuration  state  is  dominant  in  one 
region,  and  another  state  in  a  different  region,  and  so  on.  The  topological 
excitation-soliton  comes  in. 

We  now  show  that  a  soliton  solution  exists  at  %(E.  +  E  ).  For  this, 

A  D 

we  consider  the  Hamiltonian  Eq.  10-2  in  real  space.  The  potential  energy 
contributions  from  both  acceptors  and  ion  energies  can  be  lumped  into  the 
effective  interaction  and  effective  distortion  energies. 

We  easily  get  the  eigenvalue  equation  in  real  space,  corresponding 
to  Eq.  10-J: 


{^3  f  +  A(R,r)  }  4>n  =  En  <j>n 


(10-9) 


with  the  spinor  <p  =  {  U*(R,r),  V*(R,r)}  and  the  self-consistent  condition. 


A  ( R ,  r )  =  >!  (-g)l!(R,r)V*(R,r), 


(10-10) 


The  energy  difference  factor  c  may  be  written  as  e  =  (-h  /Zm^cl^d^. 
where  m^  is  the  effective  hole  mass,  R  =  (r^  +  rg)/2,  and  r  =  rA  -  rD  , 


A  B 


with  respective  hole  coordinates  r^  and  r^  .  For  simplicity,  we 
consider  a  quasi-one  dimension  in  the  direction  of  separation  of  two 


appropriate  acceptors.  Taking  the  Fourier  transformation  via  exp(ik^x^), 
where  x.  =  r,  or  R,  and  with  the  characteristic  wave  vector  k^  ,  we  get 
c  =  (  fi k  ^ / 2 )  (-  hd^)  •  Equation  10-9,  with  the  above  e,  becomes 


identical  to  the  spin  less  Dirac  equation  discussed  by  Jackiw  and  Rebbi 
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who  have  found  a  soliton  solution  for  E  =0,  with  A(x.)  =  A  tanh(x./f.) 


.mil  .1  soli  (on  Co  liavo  Fermion  numl)i*r,  equ  ivalent  ly  '^e  charge.  The  length 

parameter  (.  ran  he  eons ide red  as  a  variational  parameter.  For  sol i tons  in 
H  H  H9 

po  I  vaeet  v  I  ene  ’  ,  the  lowest  soliton  formation  energy  F,  =  (2/ti)  A  is 

S  uo 

found  when  f,  =  h  Vp/A^  .  We  now  estimate  the  soliton  formation  energy, 

Ks(r)  =  W(A  tanh  x./£  )  -  W(A)  (10-11) 

2  2  1 

with  A  =  (g  -t:  )2  from  Eq .  10-8.  Since  A (c)  <  A(0)  ,  we  expect  Eg(e) 

to  he  less  than  E  (0).  By  taking  a  mean  value  of  Eg(0)  over  the  character¬ 
istic  length,  we  easily  get  Eg(0)  =  g  ln2,  less  than  g. 

We  thus  can  have  a  soliton  solution  at  %(E.  +  E„)  ,  with  g  >  e  . 

A  U 

The  emperical  relation  of  Eq.  10-1  is  justified  in  our  physical  picture  of 
soli tons. 

We  remark  that  the  interaction  and  distortion  energies  in  Eq.  10-2 
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can  be  derived  microscopically  in  the  spirit  of  Anderson  .  When  a  distor¬ 
tion  d  is  introduced  in  the  system,  the  coupling  and  distortion  energies 

2  + 

can  be  written  by  2gAd(n  +  nD)  +  d  ,  with  n.  =  C.  C.  (i  =  A,  B) .  Minimiz- 

U  A  d  111 

ing  it  with  respect  to  d,  and  together  with  the  Coulomb  energy  U  n  n  ,  we 

C  A  d 

get  -g0(nA  +  nB)  -  g  nAng  ,  with  the  effective  coupling  g  =  2gQ  -  Uc  .  The 

first  term  can  be  lumped  into  E  and  E  without  changing  e  in  Eq.  10-2. 

A  D 

The  rest  may  be  written  after  factorizing  as 

-gnAnB  "  A  C1CB  +  A*CBCA  +  lA|2/S  (10~12> 

with  A  =  <  7.  (~g)  ctcA  >  .  The  above  is  identical  to  the  one  used  in 

D  A 

Eq.  10-2.  In  particular,  for  2g^  >  ,  the  effective  interatction  between 

n  and  n  becomes  attractive.  This  is  precisely  the  origin  of  superconduc.- 
tivitv  with  the  electron  phonon  interaction  2gn  being  greater  than  U  . 


11.0  A  GOHERENT  EXCITATION  STATE  AND  ANISOTROPIC  RAMAN 
SCATTERING  EFFICIENCIES  IN  n-  GaAs 

A  coherent  excitation  state  has  been  introduced  to  understand  the 

anisotropic  Raman  scattering  efficiencies  in  n-GaAs  under  the  resonance 

...  42 

condition 

In  the  Raman  scattering  (with  the  K  -6471A-laser)  by  coupled  plasmon 
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1.0  phonon  modes  in  n-GaAs  samples  with  n  =  3.1  -  14  x  10  cm  ,  the  Raman 
peak  due  to  the  L+-  mode  is  observed  in  the  X(Y,  Y  +  Z)  X  configuration, 

at  300  K  and  77  K,  where  X  =  (100),  Y  =  (Oil),  Z  =  (Oil),  and  X  =  (TOO) . 

However,  the  Raman  peak  due  to  the  L  -  mode  is  larger  in  the  X(Z,  Y  +  Z)X 
configuration.  With  the  K^-blbAA-laser ,  the  above  anisotropic  behavior 
is  observed  only  at  300  K. 

A  qualitative  understanding  of  the  anisotropic  Raman  efficiencies 
may  he  described  at  follows:  The  optical  scattering  efficiency  of  Raman 
processes  in  a  semiconductor,  such  as  n-GaAs,  is  known  to  depend  upon  the 
thickness,  L,  of  the  depletion  layer,  and  the  absorption  length,  &.  In 
particular,  the  coupling  of  the  plasmon  and  phonon  would  occur  only  in  the 
region  of  V.-L.  The  size  of  this  optical  active  region  would  be  sensitive 

to  the  dangling  band  and  surface  states,  since  the  band-bending  behavior 

depends  on  them. 

On  the  GaAs  (110)  surface,  the  lxl  surface  reconstruction  has  been 
observed,  and  the  2  x  2  on  the  (111)  Ga  surface,  but  only  the  1  x  1  on  the 
(111)  As  surface.  On  the  other  hand,  the  2x8  reconstruction  has  been 
reported  on  the  GaAs  (100)  surface  of  Ga ,  or  As.  Thus,  the  asymmetrical 


behavior  of  the  optical  scattering  efficiency  is  predicted  to  occur  on  the 
(100)  surface  of  GaAs,  but  not  on  the  (110)  and  (111)  planes. 


From  a  formal  point  of  viow  ,  the  Raman  scattering  tensor  in  the 

hack  scattering  X(  ,  )X  configuration  can  be  written  as 

0  0  0  D  0  0 

e 

0  0  D  +  i  0  D  0 

e 

0  I)  0  ,  0  0  D 

e 

with  1)  =  +  9u/3€  =  H  exp  (iO) 

De  =  Dp  +  °F  =  de  GXP  ’ 

where  ha/ hu  ,  ha/a  G,  ,  and  are  due  to  the  deformation  potential, 
electro-optic,  charge  density  fluctuation,  and  Frohlich  interaction 
contr Lbut ions . 

The  Raman  scattering  efficiencies  in  the  X(Y,  Y  +  Z)X  and  X(Z,  Y  +  Z)X 
configurations  become  proportional  to 

1 1)+  iDe  |  2  =  d2  +  d2  +  2dd^sin  (0-<p)  . 

in  the  non-resonance  condition,  d  »  d^  ,  we  expect  no  anisotropic 

behavior  in  the  Raman  scattering  efficiencies.  However,  under  the  resonance 

condition,  the  d  can  be  of  the  order  of  d  ,  and  a  coherent  state  comes 
e 

in;  then  the  anisotropic  behavior  in  the  Raman  scattering  efficiencies  is 
predicted  to  occur.  Tn  particular,  the  largest  anisotropy  will  occur  when 

o  -  ■}=  a/.,  . 
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We  believe  that  the  preliminary  data  can  be  understood  via  a 
coherent  state  under  the  resonance  condition.  It  is  highly  desirable  to 
carry  out  further  study  along  the  above  lines. 


12. 0  SUPERLATTICE  STRUCTURE  AND  SUPERCONDUCTIVITY  IN  HIGH 


PRESSURE  QUENCHED  CADMIUM  SULFIDE 

One  ol  tin*  intriguing  aspects  of  the  diamagnetic  anomaly  observed 

at  17  K  in  high  pressure  quenched  radinium  sulfide  is  the  possibility  of 

94  95 

having  It  i  }»li  temperature  superronduet  i  v  i  ty  ’  .  We  report  here  our 

observation  of  large  diamagnetism  in  high— pressure— quenched  cadmium 

sulfide  at  77  K,  indicating  possible  superconductivity  at  high  temperature, 
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or  a  new  condensed  state  ,  and  our  discovery  of  the  two  new  types  of 

commensurate  and  incommensurate  superlattice  structures  in  cadmium  sulfide 

via  the  X-ray  diffractometer.  97 

12.1  A  Now  Condensed  State  -  Superconductivity 

Recently,  the  evidence  for  superconductivity  at  77  K  in  high 
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pressure  quenched  CdS  has  been  reported  ’  .  The  magnetization  of 

pressure  quenched  CdS  samples  appears  to  be  like  that  of  a  type  II  super¬ 
conductor,  and  disappears  after  a  few  days.  It  is  suggested  that  this 

mi1  las  table  phenomenon  may  be  attributed  to  the  notion  of  finite  momentum 

98 

pairing,,  i.e.,  pairing  of  holes  .  We  describe  hero  our  observation  of 

such  large  diamagnetism,  not  inconsistent  with  superconductivity. 

Our  method  of  pressure  quenching  is  quite  different  from  that  of 
94,  9  5 

Homan  ot  al.  ’  whose  pressure  quenching  was  carried  out  at  a  rate 
of  10^  bar/sec,  at  room  temperature.  Our  pressure  quenching  was  carried 
out  al  77  K  bv  opening  the  pressure  oil  valve  after  the  pressure  bomb 
(in  the  range  ol  SO  -  50  Khar)  was  immersed  in  liquid  nitrogen;  then,  the 
bomb  was  disassembled.  For  the  pressure  transmitting,  fluids,  Na 0.1  and 
AgOl  were  used.  We  are  not  able  to  say  which  quenching  method  is  better. 
We  prefer  the  pressure  quenching  at  77  K,  since  a  structural  phase  transi¬ 
tion  is  known  to  occur  around  80  ~  150  K  in  CdS. 


The  start  ill}',  materials  were  pure  and  doped  (LL,  Cl,  and  other 
e  lemenl  s )  ,  CdS  grown  approximately  15  years  ago,  and  also  powder  CdS  from 
Alla  I  no  l  ean  i  e  stock  No  .  20  I  10 . 

The  measurement,  of  L  lie  ae  (11  II/.)  magnetic  susceptibility  ol  L  he 
sample  was  carried  out  by  using  pick-up  coils  and  a  lock-in  amplifier 
PAR  Model  124A.  The  reference  signal  was  picked  from  a  20  3  resistance 
connected  in  series  with  the  primary  coil.  The  inductive  part  of  the 
signal  from  the  secondary  coil  was  chosen  to  correspond  to  the  real  part 
of  the  magnetic  susceptibility  of  the  sample.  All  data  were  taken  after 
the  pick-up  coils  and  the  samples  were  immersed  in  liquid  nitrogen.  For 
the  calibration  of  susceptibility,  a  7  mg  sample  of  MnC^  was  used. 

Some  pressure-quenched  samples  were  found  to  have  a  diamagnetic 
susceptibility  of  =  (~  V.'ix)  (0.1  ~  0.2),  and  others  to  have  a  para¬ 
magnetic  susceptibility  of  y  ~  0.07.  The  magnetization  has  disappeared 
after  a  day  or  so  for  diamagnet ized  samples,  and  after  a  somewhat  longer 
time  for  pa ramagnet ized  ones. 

We  have  examined  possible  spurious  signals  due  to  the  following: 


(a)  c.oil  instability,  (b)  temperature  instability,  (c)  oxygen  in  liquid 
nitrogen,  (d)  liquid  oxygen  in  sample  holder,  and  (e)  contamination  of 
magnetic  impurities,  such  as  Fe ,  Co,  and  other  elements.  For  (a)  and 

(b) ,  the  stability  was  tested  before  taking  a  signal  from  the  samples. 

For  (c)  and  (d) ,  we  used  fresh  liquid  nitrogen,  and  tested  the  empty 
sample  holder.  (e)  is  irrelevant  for  the  diamagnetic  signal ,  and  unlikely 
to  affect  any  signals,  since  the  magnetization  has  disappeared  after  a 
few  days. 

Our  result  is  not  inconsistent  with  the  reported  superconductiv- 
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itv  in  pressure  quenched  CdS  at  77  K  ’  .  An  X-rav  examination  of  the 

samples  is  described  in  the  next  section. 


12.2  ComnuMisu rau>  and  1  iH-oniinen.su rate  Superlattice  Structures 


Cadmium  sulfide  is  known  to  have  a  wurtzite  structure  with  a  = 

4.1  56  A  and  c  =  6.713  A  at  ambient  pressure,  zinc  blende  with  a  *  5.818  A 
under  pressure,  and  rocksalt  structure  with  a  =  5.464  A  under  high 
pressure,  above  30  Kbar. 

The  low-angle  X-ray  diffraction  lines  associated  with  commensur- 


ate  super  lattice 

structures  via  an  X-ray  diffractometer 

are  observed 

1H,hkl: 

Wu  rt  z  i te 

Zinc  blende 

rocksalt 

2  d110 

2  d220 

2  dlll 

3  d110 

3  d220 

3  dlll 

2  d101 

4  d311 

2  d100 

*  d  1  01 

I'or  the  rocksalt  structure,  a  line  associated  with  the  commensurate  index 
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n  =  3  has  been  reported  .  However,  the  lines  associated  with  the  com¬ 
mensurate  indices  n  =  2  and  3,  with  hexagonal  as  well  as  cubic  structures 
are  observed  for  the  first  time,  to  our  knowledge,  in  cadmium  sulfide. 

The  X-ray  diffraction  lines  associated  with  an  incommensurate 
superlattice  structure,  via  the  X-ray  diffractometer,  are  observed  in 
high  pressure  (above  30  Kbar)  quenched  cadmium  sulfide  (black  metallic 
luster)  at 
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The  line  at  d  =  8.890  A  is  found  Lo  have  an  intensity  stronger  than  those 
of  the  regular  strong  lines  at  d  =  3.583,  3.557,  and  3.16  A  in  cadmium 
su 1 f ide . 

We  have  ruled  out,  at  the  present  time,  the  polymorphism  and 
polytypism  in  cadmium  sulfide,  since  the  strong  sharp  line  at  d  =  8.890  A 
cannot  he  present  in  those  structures.  We  conjecture  that  a  simple  fee 
structure  with  many  atoms  in  the  unit  cell  would  be  formed  in  high  pres¬ 
sure  quenched  cadmium  sulfide. 

It  is  highly  speculated  that  the  superlattice  structure  would 
make  the  effective  density  of  states  at  band  edges  increase  and  also 
make  the  number  of  degenerate  interval  leys  increase  so  that  high  tempera¬ 
ture  super— conduct ivitv  would  be  realizable 
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